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ABSTRACT 

Observations of the binary pulsar PSR J1141— 6545 using the Parkes radio 
telescope over 9.3 years show clear time-variations in pulse width, shape and 
polarization. We interpret these variations in terms of relativistic precession of 
the pulsar spin axis about the total angular momentum vector of the system 
changing our view of the emission beam. Over the nine years, the pulse width at 
the 50% level has changed by more than a factor of three, reaching a maximum 
value of nearly 13° in early 2007. Large variations have also been observed in 
the 1400-MHz mean flux density; this reached a peak of ~ 20 mJy in mid-2002 
but over the past several years has been relatively steady at ~ 3 mJy. The pulse 
polarization has been monitored since 2004 April using digital filterbank systems 
and also shows large and systematic variations in both linear and circular polar- 
ization. Position angle variations, both across the pulse profile and over the data 
span, are complex, with major differences between the central and outer parts 
of the pulse profile. We interpret the outer parts as representing the underlying 
magnetic field and fit the rotating-vector model to these regions. Modelling of 
the observed position angle variations by relativistic precession of the pulsar spin 
axis shows that the spin-orbit misalignment angle is about 110° and that the pre- 
cessional phase has passed through 180° during the course of our observations. 
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At the start of our observations, the line-of-sight impact parameter was about 
4° in magnitude and it reached a minimum very close to 0° around early 2007, 
consistent with the observed pulse width variations. We have therefore mapped 
approximately one half of the emission beam, at least out to a radius of about 
4°, showing that it is very asymmetric with respect to the magnetic axis. The 
derived precessional parameters imply that the pre-supernova star had a mass 
of about 2 M and that the supernova recoil kick velocity was relatively small, 
between 100 and 250 km s -1 , depending on the assumed systemic velocity. With 
the reversal in the rate of change of the impact parameter, we predict that over 
the next decade we will see a reversed "replay" of the variations observed in the 
past decade. 

Subject headings: pulsars: individual (PSR J1141— 6545) — relativity — radia- 
tion mechanisms: non-thermal 



Introduction 



Binary pulsars with short orbital periods and massive companions exhibit a range of 
relativistic effects, the most prominent of which are relativistic periastron precession, grav- 
itational time dilation, transverse Doppler effects, orbit decay due to loss of energy to 
gravitational waves and the Shapiro delay resulting from the passage of the pulsar sig- 
nal through the gravitational field of the c ompanion. These effects may b e described by 
the so-called "Post-Keplerian" parameters ( iDamour fc Deruelld Il985l 119861 ) and are most 
readily observed in double-neutron-star s ystems such as the Hulse- Taylor binary pulsar, 
PSR B1913+16 dWeisberg fc Tavlorl 120051 ) and the Double Pulsar, PSR J0737-3039A/B 
(IKramer et al.l 120061 ) . Another relativistic effect which is potentially observable in such sys- 
tems is the precession of t he pulsar spin axis resulting from coupling between the spin and 
orbital angular momenta (jDamour fc Ruffinilll974l : iBarker &: Q'Connelllll975l ). The angular 
rate of precession, Q p , is given by: 

q, = ^|-l 5/3 ,. mc( „ 4mp + 3mc) ,,„ (i) 



2tt 



(1 -e 2 )(m p + m c ) 4 / 3 



where G is the gravitational constant, c is the velocity of light, Pj, is the binary orbital period, 
e is the orbital eccentricity and m p and m c are the pulsar and companion masses respectively. 
For PSR B1913+16, the expected rate from Equation (CQ) is 1?21 yr _1 , corresponding to a 
precessional period of about 300 years, whereas for PSR J0737— 3039A/B the rate is about 
a factor of four larger and the precessional period is about 75 years. 

Since the emission from pulsars is believed to be beamed, probably along the open field 
lines associated with the magnetic poles on the neutron star, precession of the pulsar spin 
axis will result in changes in the beam aspect as viewed from the Earth. One would expect 
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this to lead to changes in the observed pulse profile and changes attributed to this effect 
have indee d been observed in PSR B1913+16 (jWeisberg et al.lll989l ; lKramerlll998l ) and PSR 
B1534+12 jArzoumanianlll995l : [stairs et alJl2004h . For PSR B1913+16, a significant fraction 
of the expected beam diameter has been traversed in the time since its di scovery and this has 
been used to map the two-dim ensional structure of the emission beam (jWeisberg fc Taylor 
20021 ; IClifton fc Weisbergl 120081 ) . Because of its relatively short precessional period, similar 
effects were expected to be observable fo r the Double Pulsar, but surprisingly, they h ave 
not so far been observed in the A pulsar ( Manchester et al.ll2005l ; iFerdman et al.ll2008l ). A 
possible reason for this is that A's spin-orbit misa lignment angle (the angle between th e 
spin and orbital angular momenta) is small (see also IWillems et al.ll2004[ IStairs et al.ll2006l ). 
However, large variations in the pulse profile and orbital visibility of pulsar B have been 
observed ( Burgay et al.l2005 ) and these are certainly due to some combination of relativistic 
periastron precession and precession of the B spin axis. There is now direct evidence for the 
relativistic spin precession of pulsar B. The emission from pulsar A is eclip sed for about 30 
seconds when it passes behind pulsar B (ILyne et al.ll2004j ; iKaspi et al.ll2004j) with the eclips e 
profile showing modulation related to the spin period of pulsar B ([McLaughlin et al.l 12004 ) . 
Modelling of the long-term variations in the eclipse profile are consistent with t he relativistic 



prece ssion of B's spin axis with a spin-orbit misalignment angle of about 130° (jBreton et al. 
2008h . 



PSR J1141— 6545 is a 394-ms pulsar in an eccentric 4.7-h binary orbit with a rela- 
tively massive com panion (m c ~ 1.0 M Q ) discovered in the Parkes Multibeam Pulsar Survey 
( IKaspi et al.l 120001 ). With these characteristics, relativistic effects are expected to be de- 
tectable and indeed, measurement of the relativistic advance of periastron was reported in 
the discovery paper. The expected rate of relativistic precession of the pulsar spin axis is 1?36 
yr -1 , corresponding to a precessional period of about 265 yr. Unusually for such a binary 
system, the pulsar characteristic age is relatively low, about 1.4 x 10 6 yr, suggesting that 
the initially more massive primary star evolved to form a white dwarf, the present compan- 
ion star, and in the process transferred mass to th e secondary which subsequently exploded 



leaving behind a neutr on star, the present pulsar ( Dewey &: Cordesl 



Tauris fc Sennels!l2000l ). Continued timing observations ( IBailes et al. 



9871 Kaspi et al. 2000; 



20031 : iBhat et al.ll2008f ) 



reinforced these ideas and gave detections of relativistic time dilation and orbital decay which 
were consistent with the predictions of general relativ ity. A clea r orbit al modulation in the 
timescale for interstellar scintillation was detected by lOrd et al.l (120021 ) , giving estimates of 
the orbital inclination angle i (or 180 — z) = 76?0±2?5 and the system tran sverse space veloc- 
ity, 115 ± 10 km s _1 . From measurements of the relativistic Shapiro delay, IBhat et al.l (120081 ) 
derive a value for the inclination angle of i = 73?0 ± 1?5, consistent with the scintillation 
measurement. 



Observations at frequencies around 1.4 GHz by iHotan et al.l ( I2005a[ ) showed that the 
mean pulse profile for PSR J1141— 6545 changed significantly over the five years from 1999 
July to 2004 May (MJD range 51381 - 53134), with an approximately linear increase in 
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Wio, the pulse width at 10% of the pulse peak. Furthermore, based on two measurments 
separated by about 0.7 yr, there appeared to be steepening of the gradient of polarization 
position angle (PA) dip/dcf), where ip is the PA and <p is the pulse phase (measured in the 
same units as the PA) near the profile center, at a rate of 2.3 ± 0.4 yr -1 . Variations of 
PA across pulse profiles are often well described by the "rotating-vector model" (RVM) in 
which the emission is assumed to be polarized parallel (or perpendicular) to the projected 
direction of magnetic fields in the vicinity of a magnetic pole on the rotating neutron star 



(IRadhakrishnan fc Cookdll969l ). In the RVM, ip varies as 



sin a sin [(p - <p 

tan ^ - = -r-r -r-. 7T — -r 2 

sin Q cos a — cos Q sin a cos(0 — <p ) 

ipQ is the PA at <p , the center of symmetry of the PA variation and the projected direction of 
the rotation axis for a dipole field, ( = a + /3 is the inclination of the observer's line of sight 
relative to the rotation axis of the star, a is the inclination of the magnetic axis relative to 
the rotation axis and (5 is the impact parameter of the observer's line of sight, that is, the 
minimum angle between the magnetic axis and the observer's line of sight which occurs at 
<P = 0o0 From Equation [2j 

(dip/d<p) m = sina/sin/3 (3) 
where (dip/d 4>) m is the maximum value of dip /dp at 0o- Therefore, the increase in PA gradient 



observed by iHotan et al.l ( ]2005al ) for PSR J1141— 6545 was interpreted as a decrease in \/3\, 
implying that our line of sight was moving closer to the beam axis. This and the changing 
pulse width were interpreted in terms of relativistic precession of the pulsar spin axis with a 
spin-orbit misalignment angle greater than 15° and probably less than 30°. 

In this paper we report on observations of PSR J1141— 6545 made over 9.3 years from 
1999 August to 2008 November with the Parkes radio telescope at frequencies close to 1.4 
GHz. Since 2004 April, we have observed with systems recording full polarization data 
but before that, only total intensity data were obtained. The observations and analysis 
procedures are described in Section [2] and results are presented in Section [3J In Section @] 
we describe the interpretation of the observed variations in terms of relativistic precession 
of the pulsar spin axis and the implied shape of the emission beam is described in Section |5j 
The implications of our precessional model for the formation of the system are described in 
Section EJ In Section [7] we summarize results and give our conclusions. 



1 Note that here and in Sections HI [5] and [6j ip is defined in a right-han ded coordinate system, in creasing in 
the clockwi se direction from Ea s t (loo king toward the source) following IDamour fc Tavlorl (| 1992T) (hereafter 
DT92) and Everett &: Weisberd (|200ll) . This is opposite to the IAU convention in which PA increases in the 
counterclockwise direction from North toward East looking toward the source. The polarization conventions 
are discussed in more detail in Section [21 
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2. Observations and Analysis Procedures 

We have observed PSR J1141— 6545 using the Parkes 64-m radio telescope between 
1999 August and 2008 November (MJDs 51411 to 54785) using (at different times) two 
receiving; systems, either th e center beam of the Parkes 20-cm multibeam (MB) receiver 



( IStaveley-Smith et al.lll996l ) or the "H-OH" receiver. The MB receiver has a bandwidth of 
about 300 MHz centered at about 1.4 GHz and a system equivalent flux density (S sys ) of 
approximately 29 Jy. The H-OH receiver has a wider bandwidth, approximately 600 MHz 
from 1.2 to 1.8 GHz, but a somewhat higher S sys of about 35 Jy. Both receivers have or- 
thogonal linearly-polarized feeds and a linearly-polarized, broad-band and pulsed calibration 
signal which can be injected into the feed at 45° to the two signal probes. Four different 
backend systems were used, all of which recorded data across a total bandwidth of 256 
MHz: an analogue filterbank (AFB) having 512 frequency channe ls on each polarization a nd 



recording total-intensity data using a one-bit digitiser system (see Manchester et al.ll200ll ). a 
wideband correlator (WBC) which recorded all four Stokes parameters with 1024 frequency 
channels, and two digital filterbank systems (PDFB1 and PDFB3) which also recorded full 
polarization data with 512 channels across the band. For the WBC and the PDFBs, each 
observation was preceded by a short (2-min) observation of the pulsed calibration signal. 
Observations of Hydra A, assumed to have a flux density of 43.1 Jy at 1.4 GHz, were used to 
set the flux-density scale. Off-line processing made use of the psrchive pulsar data analysis 



system (IHotan et al.l 120041 ) . 



A log of the observations is given in Table [H AFB observations, which commenced 
soon after the discovery of the pulsar, were recorded on tape and processed off-line to form 
mean total intensity profiles. Frequency channels containing known interference were given 
zero weight before data were summed across the band. For the purposes of this paper, 
observations made over intervals ranging from a few hours to a few weeks have been grouped 
as shown in Table [H Pulse widths and flux densities were determined for each observation 
and mean values and their uncertainties for the group determined from these. The flux- 
density scale for the AFB data was established by comparison with contemporaneous WBC 
observations. 

WBC and PDFB observations were typically of 1-h duration and were folded on-line 
to form 1-min sub-integrations. Frequency channels containing strong interference, affected 
by resonances in the feed or with low system gain (normally 5% of the band at the band 
edges) and sub-integrations affected by strong impulsive interference were given zero weight. 
Using the pulsed calibration signal as a reference, the variations in instrumental gain and 
phase across the band were removed, Stokes parameters formed and the data placed on a 
flux-density scale. Data were summed to form 15-min sub-integrations and mean profile 
widths and flux densities computed as for the AFB data. All profiles were recorded with 
1024 bins across the pulse period. 



Observed PAs given in this paper follow the astronomical convention, with PA mea- 
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sured from celestial North and increasing toward East (counter-clockwise looking toward the 
source). The IEEE definition of circular polarization is adopted (i.e., at a given point in space 
and looking at the source, the E- vector for a right-circular wave rotates in a counter-clockwise 
direction). Both of these definitions are in accord with IAU recommendations (Trans. IAU, 
15B, 166, 1973). Stokes V is defined in the sense left-circular minus right-circular to conform 
with established pulsar polarization conventions which pre-date the IAU recommendations. 
File header parameters defined the polarization setup for each receiver/back-end combination 
and these were used by psrchive to correct the Stokes parameters so that they conformed to 
the se definitions. For furthe r details on polarization conventions and their implementation, 
see Ivan Straten et al.l ( 120091 ). The strong millisecond pulsar, PSR J0437— 4715, was observed 



almost every session and used to check the sign of the calibrated Stoke s parameters and the 



absolute value of position angle by comparison with the results given by lNavarro et al.l (119971 ) 
which are known to be in accordance with these definitions. The linearly polarized intensity 
L = (Q 2 + U 2 ) 1 / 2 is a positive definite quantity which is biased by noise. Plo t ted va lues have 



been corrected for this bias using the relations given by lLorimer fc Kramerl (120051 ) . 



Before summing in frequency, the observed pulse profiles must be corrected for Faraday 
rotation across the band. The nominal rotation measure (RM) for PSR J1141— 6545 is — 84±2 



rad m -2 (IHan et al.l 120061 ) . RMs were measured for each epoch having polarization data by 
summing the upper and lower halves of the bandpass separately using the nominal RM, 
taking a weighted mean of the PA differences between the upper and lower bands across the 
profile, recomputing the RM and then iterating until convergence. We show below that the 
observed PA variations across the pulse profile are complex, both in time and in frequency, 
with quite different behavior in the central and outer parts of the profile. Consequently, 
RMs were determined separately for the central and outer parts; we believe that the outer 
parts represent the true interstellar RM. 

The MB receiver suffers from significant coupling between the nominally orthogonal 
feed probes. Observed polarization variations resulting from this coupling are a function 
of parallactic angle and the coupling parameters can be measured by analysis of a series of 



obser vations of a polarized source covering a wide range of parallactic angles ( Ivan Straten 



20041 ). PSR J0437— 4715 was observed every few months, typically for 10 minutes each 
hour during a 10-hour transit, for each of the WBC and PDFB configurations used with 
the MB receiver (Table H]). These data were analysed using the psrchive program PCM to 
determine the feed ellipticities and their relative orientation as a function of frequency across 
the band. Observations of PSR J1141— 6545 were then calibrated to remove the effects of the 
cross-coupling. Cross-coupling in the H-OH receiver is an order of magnitude less and this 
calibration step was not required. Figure [1] shows (truncated) PDFB1 polarization profiles 
for PSR J0437— 4715 taken with the MB and H-OH receivers using the observational and 
processing methods described above. The results are essentially identical within the noise 
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uncertainties, confirming that our calibration procedures are robust @ 



Results 



While not the primary objective of this work, it was necessary to monitor period vari- 
ations of the pulsar in order make accurate predictions of the topocentric period for use 
during t he observatio n s and to sum observations in off-line processing. Initially, the timing 
model of iBailes et al.l ( 120031 ) was used for this purpose. However, while processing the 2007 
July 18 observations, it was realised that the pulsar had suffered a sizable glitch. Glitch 
parameters result ing from a fit of a model including post-glitch exponential decay (see, e.g. 
Wang et al.l 120001 ) to pulse time-of-arrival (ToA) data from 2006 April to 2008 November 
are given in Table [2j ToAs were obtained fr om average pulse p rofiles for data segments of 
15-min duration and the program tempo2 (IHobbs et al.ll2006P l was used to determine the 
timing parameters. There was only a small relaxation of the pulse frequency toward the 
extrapolated pre-glitch solution following the glitch, i.e., the Q parameter was small. Sig- 
nificant timing noise remains after the fit and the uncertainties have been multiplied by five 
(~ yiXr)) t° allow for this. 

Table |3]gives pulsar timing parameters from a fit to post-glitch data from 2007 December 
to 2008 November; the ~ 160 d following the glitch were omitted to avoid biasing the result 
by the post-glitch relaxation. Over the fitted data span there is significant timing noise and 
the fit uncertainties have been multiplied by a factor of three. Only the pulse frequency v 
and v w ere fitted fo r ; the position and the binary parameters were held fixed at the values 
given bv iBhat et al.l (|2008l) . Parameters a re quoted in TDB-compatible units, and the "DD" 
binary model (IDamour fc Deruelldll986l ) and the Jet Propulsion Laboratory Solar-System 
ephemeris DE405 ( Standish 19981 ) were used. 



Figure [2] shows the observed time variations of the 50% mean pulse width (W^o) and 
mean flux density at 1400 MHz over the 9.3-year data span. Such dramatic long-term 
variations in pulse width (more than a factor of three) are unprecedented. The most straight- 
foward explanation is that they result from precession of the pulsar spin axis changing our 
view of the pulsar beam. The variations are complex with a ~ 400-d interval around MJD 
52800 (2003 June) where the width clearly decreases with time before increasing again. 
Recent data show that a maximum width was reached around MJD 54150 (2007 February) 
and now the width is decreasing. There is good agreement between the widths obtained with 
the different back-end systems, although the AFB widths tend to be a few percent less than 
those measured with the WBC or PDFB1 at the same or similar time owing to differences in 
the instrumental impulse response. The time variation in W\§ is similar, but with a smaller 



2 The RM for PSR J0437-4715 is only +1.5 rad m" 2 |Navarro et al.l Il997l ) and so the PA difference 
expected between f369 MHz and 1433 MHz is only about 0?4, not visible in Figure [TJ 
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relative change over the data span. 

The mean pulsed flux density has also changed dramatically over the data span with a 
broad peak around MJD 52500 (mid-2002). There is a significant day-to-day variation in the 
measured flux densities. This is indicated by the large error bars on the AFB averages (which 
typically cover a week or more) and by the scatter in the WBC/PDFB measurements (which 
typically are made on a single day), both much larger than the uncertainty in an individual 
measurement. These short-term variations can be attributed to refractive scintillation for 



which the expected timescale is about seven days (jOrd et al.ll2002l ). However, the broad rise 
and decay of the flux density over the nine-year data span is unlikely to be a scintillation 
effect and we attribute this to the changing aspect of a complex beam pattern sweeping 
across the Earth as the pulsar spin axis precesses. 

Average total-intensity profiles as a function of time are shown in Figure [3] illustrating 
the profile evolutional At early times the profile was dominated by a strong trailing peak 
with a broad ramp of emission leading up to it, whereas at late times, the profile is more 
symmetric with an approximately Gaussian shape. Figure H] shows the profile evolution in 
the form of greyscale plots with both linear (left) and logarithmic (right) intensity scales. 
The logarithmic plot clearly illustrates a striking property of the profile evolution, that the 
width at a very low level (the lowest contour is at 1% of the peak) is remarkably constant, 
suggesting that the low- level flux-density contours represent the overall beam extent. As 
will be further described in Section HJ the outer parts of the profile are relatively stable, 
both in shape and flux density. The profiles in Figures |3J and H] have therefore been aligned 
using the midpoint at a constant flux density of 2% of the maximum as a reference point. 
This is an important assumption as it affects not only the interpretation of the various 
profile components, but also the parameters derived from any timing analysis. With this 
assumption, the peak of the profile has moved in a somewhat step-wise fashion from near 
its trailing edge to near its center. This evolution will be interpreted in terms of the growth 
and decay of components representing bright regions of the beam in Section H] below. 

Figure |5] shows the mean pulse polarization profiles for PSR J1141— 6545 at two epochs 
separated by about three years. Plotted PAs refer to the band-center frequency as marked 
on the plotsjj Dramatic changes are observed over the three- year span in the (total intensity) 
pulse profile and in the polarization parameters. Once the sign differences are allowed for 
and ignoring the absolute value of PA, the 2004 profile in Figure [5] agrees well with the profile 



5 Data from some closely spaced epochs have been averaged to improve clarity. 



4 Polarization plots for PSR J1141— 6545 given by iHotan et al.l (]2005al) have the opposite sign of V and 
the opposite sense of position-angle swing to the IAU and pulsar conventions used for observational results 
in this paper. We also note that the same comments apply to th e polarization profile s for the Double Pulsar, 
PSR J0737— 3039A/B, presented bv lDemorest et al.l (|2004J) and lHotan et al.l (|2005q ). As a consequence, the 
rotation measures given in these latter papers have the wrong sign. The current best estimate of the RM 
for PSR J0737-3039A/B with the corrected sign is +112.3 ± 1.5 rad m" 2 (|Demorest et al.ll2004h . 
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given by lHotan et al.l (l2005al ) taken about two months earlier. In particular, the PA changes 
rapidly near the profile center, indicating that our line of sight traverses the beam relatively 
close to the magnetic axis. However there are clear and signficant departures from the PA 
variation expected for the simple RVM. 

The evolution of the profiles of linearly polarized intensity L and of Stokes V across 
the pulse over the 4. 5- year span of the polarization data are shown in Figure El Both show 
systematic changes as a function of time. In the case of the L profiles, components at different 
pulse phases appear to get stronger or weaker as a function of time, whereas the variation 
of Stokes V is more consistent across the whole profile. It is striking that the variation of V 
across the pulse changes smoothly from a positive to negative sign change at early epochs 
to a negative to positive change at late epochs. Especially at later times, there are clear 
dips in the linearly polarized intensity at pulse phases near ±0.01. As will be discussed 
further below, these are attributed to overlapping of approximately orthogonally polarized 
components. 

Figure [7] shows the evolution of the PA variation across the pulse. The observed PA 
variations are clearly not well described by the RVM and there are clear systematic trends in 
the PA variations as a function of time. While there are some differences between adjacent 
PA profiles which may in part result from residual calibration errors or low-level radio- 
frequency interference, in general the trends in polarization properties with time are smooth 
and consistent despite the use of different receivers and backend systems at different times, 
again showing that our calibration procedures are effective. 

Most striking is the clear difference in PA evolution between the central and the outer 
parts of the profile. Figure [8] shows the PA variations over the 4. 5- year data span averaged 
over the inner part (phase —0.0078 to 0.0068) and the outer parts (—0.0253 to —0.0117 
and 0.0136 to 0.0253) of the profile. Variations of the mean PA for the inner part of the 
pulse profile are very significant with an intial decrease followed by a rapid increase which 
slows at later epochs. For the outer parts of the profile, the mean PA at early times shows 
some fluctuations but overall is consistent with a slow decrease, flattening at later times. 
Especially at early times, the polarization in the outer parts of the profile is quite weak 
and hence the fluctuations in PA can be attributed to remaining systematic errors and/or 
contributions from the wings of the stronger emission in the center of the profile. 

These PA changes could be due to variations in RM or to aspect changes resulting from 
the precessional motion. In Figure M we show the RMs separately for the central and outer 
parts of the profile computed using the iterative method described in Section [2J For the outer 
parts of the profile there is no significant time variation. The weighted mean RM value is 
—93.4 ± 2.2 rad m -2 with a reduced \ 2 °f 3.1; the quoted uncertainty is the weighted rms 
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deviation multiplied by A^/(Xr)@ I* 1 contrast, for the central part, there is a clear systematic 
variation of measured RM with time. We believe that this is not a true variation in the 
interstellar RM - it would be much larger than RM variations in other pulsars that are 
attributed to a chang ing; path through the interstellar medium, e.g., for the Vela pulsar 



( Hamilton et al.l 119851 ) . Furthermore, it is inconsistent with the changes in the central PA 
shown in Figure [HJ Rather, we believe the observed changes in both PA and apparent RM 
for the central part of the PSR J1141— 6545 profile are due to variations in the relative 
amplitudes of observed profile components with different PAs and spectral indices resulting 
from the changing aspect of the line of sight relative to the beam axis as the pulsar spin 
axis precesses. This is somewhat analo gous to the apparent RM var iations (as a function of 



pulse phase) seen in PSR B2016+28 by iRamachandran et al.l (120041 ) which are attributed to 
over-lapping non-orthogonal components. 

We adopt PAs in the outer zones as representing the underlying magnetic field structure. 
Consequently, these PA variations should be consistent with the RVM. The two sides are 
consistently at approximately the same PA, as expected in the RVM. Also, in some cases, 
e.g., the 2007 August 4 profile shown in Figure El there is evidence for the PA variations 
expected for the RVM at the inner edges of the outer PAs. Since there is no signficant RM 
variation for these parts, we interpret the observed slow decrease in outer-zone PAs shown 
in the right panel of Figure |8] as intrinsic and resulting from precession of the pulsar spin 
axis. 



4. Modelling of the Precessional Changes 

By modelling the observed PA variations in terms of precessional motion of the pulsar's 
spin axis, we can infer additial properties of the system, for example, the spin-orbit misalign- 
ment angle and the precessional phase. As a first step in the modelling, we fit the observed 
total-intensity profiles (Figure E]) with gaussian components. A total of six components are 
fitted, although for later epochs only four are required. We have used fixed central phases 
for the components; only the amplitude and width of each component is allowed to vary as 
a function of epoch. The component phases, given in Table HJ were chosen using an iterative 
process by fitting the first and last observed profiles. Fitted components and their sum are 
shown in Figure [TD] for a selected set of profiles which cover most of the observed data span, 
showing that an excellent fit to the observed profiles is obtained. This retroactively justifies 
our decision to fix the component phases. Figure [11] gives the widths and mean flux densi- 
ties of the fitted components. The mean flux densities are scaled so that their sum equals 
the smoothed flux-density variation shown in Figure [2j For components 3 and 5, the fitted 
amplitudes dropped to zero after MJDs 52100 and 52500, respectively. 



5 This RM value is used to rotate observed PAs to the required reference frequency. 
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lt is striking that the amplitude of each component varies smoothly from one epoch to 
the next and that the widths of the components are relatively stable. This suggests that 
these components represent real physical zones of emission which are fixed in longitude on 
the star. We note that most of the observed flux density variation comes from the central and 
trailing components. Although the leading component 1 is relatively weak at early epochs 
and dominates the profile at late epochs, its absolute amplitude remains approximately 
constant across the whole data span. 

As discussed above, the observed polarization PA and L variations (Figs [7] and [6]) 
strongly suggest that these components may be grouped into two overlapping and quasi- 
independent emission zones, one dominant in the wings of the observed pulse and the other 
dominant in the central regions. We note that similar non-RVM variations (not accountable 
for by overlapping orthogo nal polarization modes) are seen in other shor t-period pulsars, for 



example, PSR B1913+1 6 (IBlaskiewicz et al.lll99lf ) and PSR B1534+12 flArzoumanian et al. 



19961 : IStairs et al.ll2004j ). The outer components (1 and 6) may be identified with the outer 
zones in the PA and L plots. As discussed above, the PAs in these zones together are well 
fitted by a single RVM. Note that, as discussed in Sections [H and [2} for modelling purposes, 
we have adopted a right-hand coordinate system in which the sign of PA is reversed com- 
pared to the astronomical convention used in Section [31 We also choose O = 0.0 and hold 
this fixed for all epochs. 

Most previous discussions of relativistic spin precession in this and other pulsars (e.g., 
Weisberg et al.lll989t lKramerlll998t IStairs et al.ll2004[ iHotan et al.ll2005af ) have concentrated 
on measurement and interpretation of the rate of change of impact parameter. For PSR 
J1141— 6545, the impact parameter is not very well determined because of the presence of 
the evidently independent emission in the central part of the profile. However, since we have 
measured absolute position angles, we have another, potentially more sensit ive, observable: 
the central PA of the RVM, ip (cf. Istairs et al.ll2004[ iKramer k Wexll2009l ). In the RVM, 
this is the direction of the spin axis of the star projected on the sky plane (see Figure [T2|) . 
This of course varies as the pulsar precesses: 



IpO = ^asc + V 



(4) 



where Q asc is the longitude of the ascending node and rj, the precessional longitude, is the 
angle of the pulsar spin axis relative to the ascending node projected on the sky plane. In 
this equation, ipo is the observed PA corrected to infinite frequency using the measured RM 
of —93.4 rad m -2 . The polar angle A varies with the precession according to: 

cos A = cos 5 cos i — sin 6 sin i cos $ (5) 

cos rj sin A = sin 5 sin $ (6) 

where 6 is the spin-orbit misalignment angle, $ is the precessional angle measured in the 
orbit plane from — j and $o is the value of $ at t — to (DT92) and 

$ = fi p (t-t ) + $ . (7) 
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From Equation (j3J), since f2 asc is effectively constant^, the precessional variation in ipo is 
determined by the variation of rj. From Equation (jSJ), we can deter mine cos A a s a fu nction 



of time and the fixed angles 5, i and $o- We take i = 73?0 from iBhat et al.l (120081 ). take 
t — to at MJD 53000.0, near the center of the data span and the beginning of our polarization 
monitoring, and take the GR prediction for the spin precession rate fl p , approximately 1?36 
yr _1 . Then, using Equation [6] (where there is no sign ambiguity in sin A since < A < 7r), 
we obtain cos 77. From the same spherical triangle used to derive Equation we also have 

cos 5 = cos A cos i — sin i sin A sin 7] (8) 

giving us sin rj and hence rj as a functio n of time (and the fix ed but unknown angles Q aS c, 8 



and $ ) without further ambiguity (cf. Kramer fc Wexl 120091 ) . 



Consequently, the set of three parameters, 5, $ , A^o> (where Aip absorbs the con- 
tribution to ipo by f2 asc and any other intrinsic constant offset to the PA values) describes 
the variation of tpo as a function of time. At the same time, the PA swing at each epoch is 
ideally described by the RVM, which uses two further parameters: the magnetic inclination 
a and the impact parameter /3. (We have held the central longitude of the RVM, </> , at zero 
as discussed above.) However, since 

A = 7r-C = vr-a-/3 (9) 

the angle /3 can be computed from Eqn. (jSJ) for a given a. As a result, only four parameters 
should describe the behavior of the PA variations at all epochs. 

As described above, the central range of the PA swing shows clear deviations from an 
RVM. However, in a blind search for the best fit solution, it nevertheless adds valuable 
information as it clearly indicates (by its slope) that /3 < (in the DT92 convention). It can 
also be noted that the slope steepens and then flattens again towards the end of our dataset 
(see Figure [7|). In order to use the geometric information provided by the central part but 
to still tie the results to the outer wings, we decided to decrease the relative weight of the 
central part in the fitting process (by multiplying its errors by a factor of five@). We allowed 
a PA offset between the central part relative to the RVM determined by the outer wings, 
and this was fitted for independently for each epoch. We take longitudes in the range of 
|0| < 3° (0.00833 in phase) to define the central region, whereas PAs in the transition zone 
3° < |0| < 5° (0.00833 to 0.0139 in phase) were ignored in the fit. 

As a first step of the fitting procedure, we constrain the allowable values of S and $0 
using just the outer- zone PAs. For this purpose, we compute the mean PA of the outer wings 



6 The orbit does precess since it is the total angular momentum which is conserved. However, since the 
orbital angular momentum is much greater than the spin angular momentum of either star, the variation is 
very small and can be neglected in this analysis. 

Varying this factor between three and ten did not have a major impact on the results described below. 
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and fit its behavior as a function of time (see Figure [T3|) . We use a grid search in the 5 - 
$o plane, where at each grid point we compute the y 2 for the fit of the model values (based 
on Equations S] - E]) to the computed values of -0o (Figure [TBI) . Figure IT41 is a plot of x 2 i n 
the 5 - $o plane, showing two clear and well defined solutions. For this plot, we have taken 
i = 73?0; the plot for 180 — i is mirror-symmetric. 

The solution with $ ~ corresponds to A > tt/2 while the other solution implies 
A < 7r/2. We can use the information provided by the PA swing to identify the appropriate 
solution in the following way. From the PA swing (in the DT92 convention), the slope is 
negative and hence < 0. Also, the fact that the PA becomes flat at the outer wings suggests 
an outer line of sight (that is, the line of sight is on the equatorial side of the magnetic axis 
when <f) = 0o, see, e.g., Lorimer & Kramer 2005) rather than an inner one. The combination 
of both is only possible in a solution with £ = ct + (3 > 7r/2, corresponding to solution where 
A = 7r — ( < tt/2. This uniquely identifies the region around 5 ~ 90° and $o ~ 175° as the 
correct solution. 

The best solution from this fitting procedure corresponds to A being very close to zero. 
We also know that (3 is close to zero from the rapid swing of PA in the central region and the 
flat PAs in the outer zones which are only a few degrees from the pulse center. Therefore, 
from Equation (Q, we have that a ~ n. However, the total-intensity profiles and PA swing 
show that the magnetic axis cannot be exactly aligned (or counter-aligned) with the rotation 
axis. The fact that the sign of the central PA swing remains the same at all epochs also 
shows that /3 does not change sign as a result of the precessional motion. This suggests 
that the fit to the PA data as a whole will lie somewhat to the right of the best-fit point 
shown in Figure HH This idea is confirmed by a blind fit to all PA data using only the four 
parameters 5, Aip and a. The search is again performed on a S - $o grid, where at each 
grid point we execute a simplex algorithm to minimize for the other two parameters. The 
results (collapsed in Aip and a space) are shown in Figure [Tj)J We show in grayscale the 
previous results of the PA-offset fit and overlay contour levels for the global RVM fit to the 
PA data. The model variation of ipo, shown in Figure [131 has a reduced \ 2 of 7.4 with 920 
degrees of freedom. Assumed and derived parameters for the system are given in Table [5j 
Figure [TBI shows the joint probability density function for the spin-orbit misalignment angle 
S derived from the fits to the PA variations and the global RVM fit (Figure [TBI) . 

Given the geometry of the system we can obtain the model time variation of impact 
parameter (3 using Equations [5] - [7] and [9j The result, given in Figure [TTl shows that the 
observer's line of sight has approached the magnetic axis direction through most of our 
observed data span, reaching its closest approach around MJD 54000, and is now receding 
again. This reversal in the variation of f3 is supported by the fact that the pulse width 
reached a maximum at about MJD 54200 and is now decreasing (Figure [2]). Note that we 
have obtained these results by fitting only the polarization data, i.e., we did not make use of 
the total intensity information. Although the over-all variation of pulse width is consistent 
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with the precessional model, the detailed variations are more complex. This shows that the 
assumption of a uniform circular beam is not appropriate. 



5. Shape of the emission beam 



Given the variation of /3 over the data span (Figure [T7I) and the time variation of 
component shapes and amplitudes (Figure fTTl) . we can compute the two-dimensional intensity 
profile of the emission beam over the traversed region. For each epoch we compute the 
traverse of our line of sight across the polar region and accumulate the pulse intensity in a 
two dimensional grid centered on the magnetic axis. The resulting beam pattern is shown in 
Figure [181 Because of the reversal in the time-derivative of (3 when it was close to zero, we 
only see one half of the polar-cap region. Despite this we can clearly see that the beam is 
quite asymmetric with no evidence for a core-cone or ring structure that is symmetric about 
the magnetic pole. The partially filled beam can be described as "patchy", albeit with just 
one major patch in the region scanned so far. Although this is the first two-dimensional map 
of an emission beam to clearly show such patchy structure, there is good evidence that the 



emission beam from mos t pulsars is best described in this way ( ILyne fc Manchester! 11988 
Han fc Manchester! EoOlh . 



Although the observed pulse width is about average (the median pulse width for all 
pulsars is about 10°), the small value of 180° — a ~ 20° implies that the intrinsic beamwidth 
is small. Figure [18] shows that the emitting region fits within a circle of radius about 4° 
centered on the magnetic axis. Excluding millisecond pulsars, observed pulse widths are 
generally consistent with the relation 



P 



6?5P 



-0.5 



(10) 



flGil et al.lll993l : iKramer et al.lll994f ). For PSR J1141-6545 the predicted value of p is about 
10?3, much larger than the emitting zone traversed so far. 



PSR J1141— 6545 was not detected in the Parkes 70cm survey ([Manchester et al.lll996l ) 
although, even at the present relatively low flux-density levels, a detection with signal/noise 
ratio of the order of 50 would have been expected. Observations within half a beamwidth 
of the pulsar position were made 1993 July 14 (MJD 49182). Figure [17] shows that the 
impact parameter (3 at that time was about —8°. This non-detection therefore suggests 
that the beam half- width in latitude is < 8°, although this must be qualified because of 



the patchy beam structure. There are good arguments (e.g., iNarayan fc Vivekanandlll983 
Manchester! Il996l ) that beams in young pulsars are elongated in the latitude direction and 
hence more fan-like than circular. An elongated beam has a lso been suggested for PSR 



B1913+16 flWeisberg & Tavlorll2002l ; IClifton fc Weisberd 120081 ). Observations over the next 
decade or two will establish whether or not this is the case for PSR J1141— 6545 as /3 returns 
to large (negative) values. 
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6. Implications for the progenitor system 



There is good and increasing evidence (e.g. iJohnston et al.ll2005l ; iNg &: Romanil 120081 ) 
that the velocity and spin vectors of young pulsars are nearly aligned. This suggests that 
the kick from the supernova (SN) explosion, the dominant f actor in determining; the pulsar 
space motion, may also determine the spin direction (see also lSpruit fc Phinneylll998l ). Since 
PSR J1141— 6545 is a young pulsar, its spin direction is therefore likely to be aligned with the 
SN kick. While precession will have caused the azimuthal angle of the pulsar spin relative to 
the (post-SN) orbital angular momentum to wrap many times over the lifetime of the pulsar, 
the polar angle (5 in this paper) should still reflect the direction between the kick and the 
post-SN orbit normal. 

We have explored if this expectation can set any interesting constraints on the kick 



(1996 


) and 


Wex et al. 


(2000) 



cos 5 



We follow the pro cedure set out in iThorsett et al.l (20051) a n d the Bayesian version in 
Stairs et al.l (120061 ) . which follows the principles of I Wex et al.l ( 120001 ) but which puts pri- 
ors on the system velocity, the angle 6 between the pre- and post-SN orbit normals^ and 
the orientation Q of the binary system on the sky. Each sampling of the priors determines 
the coefficients of a quadratic equation in m 2 j, the mass of the pre-SN star (Equation 8 of 



Thorsett et al.ll2005l ). which may or may not have a solution that falls within acceptable 
ranges for the pre-SN mass and the size of the (assumed circular) pre-SN orbit. As this is 
a young system, we did not evol ye the system's motion back in time through the Galaxy 
to identify plausible birth sites ( [Wex et al.l 120001 ). nor did we evoly e the orbital size and 
eccentricity to account for gravitational radiation losses (JPetersJ |l964J) , but simply assumed 
the current velocity and orbital parameters reflect the birth properties. Note that the re- 
sults are therefo re independent of the assumed distance; we assumed a distance of 3.7 kpc 
( lOrd et al.l 120021 ). We put uniform priors on 9 (sampling from 0° to 180° and testing both 
positive and negative values in each trial) and Q (0° to 360°). We sampled the companion 
mass over (1.02±0.01) M Q and th e total system mas s over (2.28911 ±0.00026) M , based on 
the orbital parameters reported in lBhat et al.l (120081 ) . We subtracted these to get the pulsar 
mass for th e given trial, and used the mass function computed from the orbit presented in 
Bhat et al.l ( 120081 ). 0.176550265 M Q , to derive sinz. We considered only the cases where 
cos i > in order to match the modeling in the previous section. 



While there is a scintillation velocity measurement for this system ( lOrd et al.l l2002h 
which could in principle provide some constraints on the angle between the proper motion 



8 Note that this is the angle labeled 5 in IThorsett et al.l (|2005l) and IStairs et al.l (|2006l ) and corresponds 
to the spin-orbit misalignment angle for a recycled pulsar whose companion has undergone a SN explosion. 
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and f2, similar measurements have been shown to be quite unreliable i n the Double Pulsar, 
probably because of the effects of anis otropics in the interstellar medium (IRansom et al.ll2004l ; 
Coles et al.ll2005l ; Kramer et al.ll2006l ) and hence we do not make use of this information here. 
Instead we investigate Maxwellian velocity distributions with 1-dimensional dispersions of 
50 km s _1 and 100 km s _1 (these dispersions are of order the transverse velocity derived from 
the scintillation measurements), combined with a second Maxwellian with 1-dimensional 
dispersion of 12 km s _1 representing the pre-SN peculiar velocity. We then deri ve Bayesian 



proba bility distributions for the pre-SN and kick parameters as described in IStairs et al. 



( 120061 ). considering pre-SN masses in the rang e 1.4 - 10.0 M to allow for the binding energy 
of the neutron star (jLattimer fc Prakashll200ll ) . For each acceptable solution of t he quadratic 



equat ion, the relevant set of parameters is assigned a uniform likelihood (as in IStairs et al. 



20061 ) for a "natural" weighting. We then make use of the probability distribution for the 
post-SN pulsar spin-orbit misalignment angle 5 obtained from the PA fitting (Figure ITBj) to 
give a constrained weighting on the range of acceptable pre-SN parameters and kick velocity. 
The derived probability density functions for the pre-SN stellar mass, the kick velocity and 
the spin-orbit misalignment angle after constraining with the allowed misalignment-angle 
distribution (Figure [To 7 ]) are shown in Figure dHJ Table gives the median values of the 
parameters and the limits at 68% and 95% confidence. We note that the natural weighting 
tends to disfavour spin-orbit misalignment angles around 90°, so the constrained distribution 
is biased toward slightly larger angles compared to Figure | 



Not unexpectedly, the assumed distribution of post-SN system velocities has a significant 
effect on the likely kick velocity, with larger values requiring a larger kick. However, the most 
probable kick velocity is relatively small in both cases, about 100 km s _1 for a dispersion of 
50 km s" 1 and about 180 km s _1 for a dispersion of 100 km s _1 . Median velocities (Table [6]) 
are somewhat higher since the distributions have a high- velocity tail. Never-the-less, they are 
still small co mpared with mean or median of the inferred three-dimensional pulsar vel ocity 
distribution ( jLyne fc Lorimerlll994l ; iHobbs et al.ll2005l ; iFaucher-Giguere fc Kaspil 120061 ) . 



Relatively low progenitor masses are strongly favoured with median values of about 
2 Mm- The puls ar mass implied by the timing fits f or companion mass and total mass 



( IBhat et al.l 120081 ) is just 1.27 M , lower than average ( IThorsett fc Chakrabartylll999l ). but 
very similar to the mass of PSR J0737— 3039B, the second-born and slow pulsar in the 
Double Pulsar system. This system also has a relatively low implied kick velocity from 
the B pulsar formation (e.g . Piran fc Shavivl l2005t IWillems et al.l l2006t IStairs et al.l 120061 ) 



but see also iKalogera et al.l ( 120081 ). These low masses and kick velocities may imply for- 



dwarf ( 


Podsiadlowski et al. 


2005 


) but other evolutionary histories are also possible 


Dewi & Pols 


2003 


)■ 



[e.g., 



The inferred spin-orbit misalignment angle 5 in the PSR J1141— 6545 system is large, 
with a most probable value of between 100° and 110°. That 5 is not small (in contrast to 
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pulsar A of the Double Pulsar system, IStairs et al.ll2006l ; iFerdmau et al.ll2008f ). is consistent 
with the major changes that we see in the pulse profile as a result of the precession and in 
fact results from the fitting of the precessional model to these changes. 



7. Summary and Conclusions 

The dramatic long-term changes in pulse shape and amplitude observed for PSR J1141— 6545 
(Figures [2] and [3]) are unpreced ented in pulsar as tronomy. Pulse profile changes are observed 



in mode- changing p ulsars (e.g. IWang et al.ll2007l ) and in radio emission from magnetars (e.g. 



Camilo et al.ll2007l ) but these are of a quite different character to the variations seen in PSR 
J1141— 6545 and are certainly caused by fluctuations in the emitted pulsar beam. Slow and 
systematic variations similar to those observed in PSR J1141— 6545 are seen in the pulse 
profiles of PSR B1913+16 and PSR B1534+12, but these are of much smaller amplitude. 
All three of these pulsars are in close binary orbits with massive companions. Consequently, 
the observed profile variations most probably result from precession of the pulsar spin axis 
resulting from spin-orbit coupling changing our view of the pulsar beam. 

Precession of the pulsar's spin axis has resulted in large and systematic changes in both 
the amplitude and the shape of the observed pulse profile as our line of sight scanned across 
the polar region. There are clearly large variations in emissivity across the nominal beam 
area, consistent with the idea of a partially filled or patchy beam. The large variations in 
the observed mean pulse profile raise the issue of profile alignment when comparing different 
epochs. Figure 0] shows that the pulse width is relatively stable at the lowest contour levels. 
We have chosen to align the profiles according to the mid-point of a constant flux-density level 
near their extreme wings. This choice is motivated by the fact that the outer components 
are relatively stable, both in flux density and shape (Figure [TT1) . With this choice, alignment 
of the polarization features over the data span is also very stable, giving added support to 
our method of defining the absolute pulse phase relative to the star. We note that this 
choice also defines a reference phase for pulse timing. It is interesting to note that, for PSR 
B1913+16 also, the pulse width defin ed by the lower beam con tours is more constant than, 



say, the pulse width at the 50% level (IWeisberg fc Taylorl 120051 ). 



Major changes are also observed in the pulse polarization parameters over the 4.5 years 
in which the polarization has been monitored. There has been a steady evolution of the 
circular polarization with the sign of the sense reversal near the profile center changing 
from positive-to-negative to negative-to-positive over the data span. Variations in the linear 
polarization are complex with quite different behavior in the central and outer zones of the 
pulse profile (Figure [7]). Our choice of taking the outer- zone PAs to represent the underlying 
magnetic-field structure is supported by a number of factors. First, the outer-zone PAs are 
well fitted by the rotating vector model (RVM) whereas this is less true of the central-zone 
PAs. More importantly, there is a large and variable offset between the central PAs and the 
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outer-zone PAs. The time variation of the mean central-zone PA is rapid and inconsistent 
with precessional motion. Furthermore there is an apparent time variation of RM for the 
central zone which is not present for the outer zones. We attribute the complex behavior 
of the central-zone polarization to the overlapping of independent components which have 
different PAs and different spectral indices. It may be possible to model the polarization 
variations based on the variations of the total- intensity components shown in Figure [11] but 
this is beyond the scope of the present work. We note that similar but much less dramatic 
non-RVM variations are observed in the central parts of the profiles for PSR B1913+16 and 
PSR B1534+12. 

Fitting of a precessional model to the observed PA variations leads to the conclusion 
that the spin-orbit misalignment angle 5 is very large in this system, with the pulsar spin 
axis nearly orthogonal to the orbital angular momentum vector. This is consistent with 
the large profile shape and polarization changes as the spin axis precesses. Unfortunately, 
because of the complex and poorly understood emission physics, it is not possible to turn 
the argument around and use these observations as a test of the precessional predictions of 
general relativity. It is interesting to note though that a large spin-orbit misalignment angle 
is quite possible despite the low pro genitor mass and m odest kick velocity. This result is 



not in conflict with the conclusion of iHotan et al.l (l2005al ). that the spin-orbit misalignment 



angle is probably less than 30°, since that was on the proviso that the pulsar was not at a 
precessional phase when the angle between the observer's line of sight and the spin axis (A 
in our terminology) was changing rapidly. In our solution, A is close to zero and is changing 
rapidly. 

Our observational data span encompassed the time when the precessional longitude 
passed through 180°. This means that the impact parameter (3 reached an extremum during 
our data span. In fact, reached a minimum value very close to zero at about MJD 54000 
(Figure [TTj) . This solution is consistent with the observed maximum in the pulse width at 
around the same time and implies that we will only ever traverse one side of the emission 
beam. The slope of the (3 time-variation has now reversed and so we are now retracing our 
earlier path across the beam. We therefore predict that over the coming decade we will see 
a time-reversed "replay" of the recent profile amplitude and shape evolution. 
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Table 1. Observations of PSR Jl 141 -6545 



Mean 


Mean 


MJD Range 


Receiver 


Backend 


Ctr. Freq. 


Nr of 


Int. Time 


Date 


MJD 






System 


(MHz) 


Observations 


(h) 


1999/08/21 


51411.7 


51407-51413 


MB 


AFB 


1390 


11 


1.72 


1999/10/02 


51452.9 


51451-51454 


MB 


AFB 


1390 


8 


1.23 


1999/12/01 


51513.5 


51498-51529 


MB 


AFB 


1390 


4 


0.66 


2000/02/04 


51577.7 


51554-51632 


MB 


AFB 


1390 


12 


2.10 


2000/06/27 


51721.6 


51710-51754 


MB 


AFB 


1390 


4 


0.99 


2000/11/23 


51851.6 


51841-51940 


MB 


AFB 


1390 


5 


0.80 


2001/04/24 


52022.9 


51969-52115 


MB 


AFB 


1390 


8 


0.78 


2001/10/20 


52201.9 


52132-52251 


MB 


AFB 


1390 


6 


0.60 


2002/03/28 


52361.0 


52305-52428 


MB 


AFB 


1390 


6 


0.60 


2002/07/30 


52485.0 


52459-52507 


MB 


AFB 


1390 


8 


1.14 


2002/11/21 


52598.9 


52571-52624 


MB 


AFB 


1390 


6 


0.60 


2003/03/17 


52715.4 


52660-52770 


MB 


AFB 


1390 


4 


0.40 


2003/11/19 


52961.6 


52914-53004 


MB/H-OH 


AFB 


1390 


6 


0.52 


2004/04/07 


53102.5 


53102-53102 


H-OH 


WBC 


1375 


1 


1.00 


2004/04/23 


53118.0 


53103-53146 


H-OH 


AFB 


1390 


3 


0.30 


2004/07/04 


53190.4 


53183-53194 


H-OH 


AFB 


1390 


5 


3.20 


2004/07/07 


53193.3 


53193-53194 


H-OH 


WBC 


1375 


8 


8.00 


2004/08/31 


53248.2 


53223-53282 


H-OH/MB 


AFB 


1390 


3 


0.30 


2004/10/31 


53309.6 


53309-53311 


MB 


WBC 


1433 


12 


11.76 


2004/11/01 


53310.2 


53306-53311 


MB 


AFB 


1390 


8 


3.85 


2005/01/02 


53372.1 


53371-53372 


MB 


WBC 


1433 


3 


3.00 


2005/01/02 


53372.5 


53371-53372 


MB 


AFB 


1390 


6 


5.70 


2005/04/23 


53482.9 


53482-53484 


MB 


AFB 


1390 


9 


8.55 


2005/06/05 


53526.1 


53526-53526 


MB 


WBC 


1433 


2 


2.00 


2005/08/04 


53586.1 


53586-53586 


MB 


WBC 


1433 


12 


7.70 


2005/08/13 


53594.6 


53522-53620 


MB 


AFB 


1390 


4 


2.16 


2005/09/05 


53618.7 


53617-53620 


MB 


PDFB1 


1433 


3 


2.05 


2005/12/27 


53730.8 


53714-53741 


MB 


AFB 


1390 


3 


1.55 


2006/01/05 


53740.8 


53735-53745 


MB 


PDFB1 


1433 


2 


1.95 


2006/02/04 


53770.5 


53770-53770 


MB 


PDFB1 


1433 


2 


1.52 


2006/03/03 


53797.7 


53797-53797 


MB 


PDFB1 


1369 


2 


1.67 


2006/04/09 


53834.5 


53834-53834 


MB 


PDFB1 


1369 


2 


1.67 


2006/05/10 


53865.5 


53865-53865 


MB 


PDFB1 


1369 


2 


1.67 


2006/05/28 


53883.5 


53883-53883 


MB 


PDFB1 


1369 


2 


1.67 


2006/07/10 


53926.7 


53921-53932 


MB 


PDFB1 


1369 


2 


1.37 


2006/08/14 


53961.3 


53961-53961 


MB 


PDFB1 


1369 


2 


1.67 


2006/09/27 


54005.5 


53995-54015 


MB 


PDFB1 


1369 


4 


2.83 


2007/01/19 


54119.8 


54119-54120 


H-OH 


PDFB1 


1433 


4 


3.00 


2007/02/09 


54140.7 


54137-54147 


H-OH 


PDFB1 


1433 


3 


1.63 


2007/05/03 


54223.5 


54223-54223 


H-OH 


PDFB1 


1433 


4 


3.33 


2007/06/06 


54257.4 


54257-54257 


MB 


PDFB1 


1369 


2 


1.67 


2007/07/18 


54299.2 


54299-54299 


MB 


PDFB1 


1369 


2 


2.13 


2007/08/04 


54316.0 


54315-54316 


MB 


PDFB1 


1369 


2 


2.13 


2007/11/05 


54409.1 


54409-54409 


MB 


PDFB1 


1369 


2 


2.13 


2008/01/24 


54489.6 


54489-54489 


MB 


PDFB1 


1369 


1 


1.06 


2008/02/24 


54509.6 


54509-54509 


MB 


PDFB1 


1369 


1 


1.06 


2008/08/03 


54681.2 


54681-54681 


MB 


PDFB3 


1369 


2 


2.11 


2008/11/14 


54784.8 


54784-54784 


MB 


PDFB3 


1369 


2 


2.13 



- 25 - 



-26 - 



Table 2. PSR J1141— 6545 glitch parameters 



Parameter 


Value 


Glitch epoch (MJD) 


54277 ± 20 


Au g /u 


(5.890 ±0.006) x 10- 7 


AVgjV 


(5.0 ±0.9) x 10~ 3 


Q 


0.0040 ± 0.0007 


r d (d) 


495 ± 140 


Data span (MJD) 


53834 - 54785 


Number of ToAs 


254 


Rms residual (/xs) 


220 



Table 3. PSR Jl 141— 6545 post-glitch timing parameters 



Parameter Value' 



R.A. (J2000) 


n h 41 m Q7:0140 b 


Dec. (J2000) 


-65° 45' 19'/1131 b 


Pulse Frequency {v) (Hz) 


2.538723048486(4) 


Pulse Frequency time-deriv. (z>) (s~ 2 ) 


-2.77615(8) x 10~ 14 


Epoch (MJD) 


54637.00 


Dispersion Measure (cm -3 pc) 


116.080 b 


Binary Period (d) 


0.1976509593 b 


Binary Period time-derivative 


-4.3 x 10~ 13b 


Orbit semi-major axis (s) 


1.858922 b 


Eccentricity 


0.171884 b 


Periastron time (MJD) 


51369.8545515 b 


Longitude of periastron (°) 


42.4561 b 


Longitude time- derivative (° yr _1 ) 


5.3096 b 


Grav. redshift /time dilation (7) (s) 


0.000773 b 


Data span (MJD) 


54435 - 54785 


Rms timing residual (//s) 


56 



Uncertainties in the last quoted digit are given in parentheses. 
b From iBhat et al.l J2008f ) . 
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Table 4. Central phases for the profile gaussian components 



Component 1 2 3 4 5 6 



Central phase -0.00892 0.00264 0.00642 0.00822 0.01325 0.01667 



Table 5. Model parameters from the RVM fit to the observed posit ion- angle variations 



Parameter Value a 



Assumed parameters: 

Orbit inclination % 73° 

Precession rate il p 1?36 yr -1 

RVM reference phase </> 0.0 

Reference time for precession MJD 53000.0 
Derived parameters: 

Spin-orbit misalignment angle 5 93° (—9°, +16°) 

Precession angle $ 175?5 (-2?8,+l?8; 

Offset in PA A^ — 118?7 ± 0?2 

Magnetic inclination a 160° (-16°, +8°) 



Most probable value and 68% confidence limits. 



Table 6. Median values and confidence limits for parameters derived from the constrained 

model distributions 



Parameter 


Vel. Disp. (km s x ) 


Median 


68% 


limits 


95% limits 


Progenitor mass (M Q ) 


50 


1.80 


1.54 


- 2.10 


1.42 - 2.48 




100 


2.24 


1.58 


- 3.06 


1.42 - 4.00 


Kick velocity (km s _1 ) 


50 


116 


86 


- 164 


56 - 240 




100 


258 


150 


-434 


94 - 452 


Spin-orbit misalignment angle (°) 


50 


116 


101 


- 131 


92 - 149 




100 


106 


94 


- 120 


87 - 135 
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Pulse Phase Pulse Phase 

Fig. 1. — Mean pulse profiles and polarization parameters for PSR J0437— 4715 at two 
epochs, 2006 October 03 (MJD 54011) using the multibeam receiver (left) and 2007 January 
18 (MJD 54118) using the H-OH receiver (right). The respective integration times were 
3.19 h and 2.13 h and the observing and data analysis methods were identical to those used 
for observations of PSR J1141— 6545. In the lower panel of each plot the solid line is the 
total intensity (Stokes J), the dashed line is the linearly polarized intensity and the dotted 
line is Stokes V. To better illustrate the details of the polarization variation, the plots have 
been truncated at 0.15 of the peak I amplitude. The position angle of the linearly polarized 
component is shown in the upper panel of each plot. Note that postion angles are defined 
according to the IAU convention; they are absolute and apply to the frequency marked on 
each plot. 
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Fig. 2. — Mean 50% pulse widths and mean flux densities at 1400 MHz versus time for PSR 
J1141— 6545. AFB data are marked with crosses and WBC and PDFB1 data are marked 
with dots. Error bars are ±1 standard deviation of the mean. The curve is a spline fit to 
average flux densities which represents the probable intrinsic variation. 
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Fig. 3. — Mean total intensity (Stokes /) pulse profiles for PSR J1141— 6545 normalised to 
the same peak amplitude. The mean MJD for each profile is shown (note that time increases 
downward). The procedure for phase alignment of profiles is described in the text. 
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Fig. 4. — Greyscale plots showing the time variation of the mean pulse profile for PSR 
J1141— 6545. The greyscale linearly interpolates between the observed profiles. For the left- 
hand plot the profiles have been normalised to unity peak amplitude (as for Figure [3]) and 
the greyscale is linear with contour lines at 0.1, 0.3, 0.5, 0.7 and 0.9 of the profile peak. 
For the right-hand plot the profile amplitudes are in flux density units and the greyscale is 
logarithmic. The contour lines are at 0.01, 0.02, 0.05, 0.1, 0.2 and 0.5 of the largest profile 
peak. 
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-0.05 0.05 -0.05 0.05 

Pulse Phase Pulse Phase 

Fig. 5. — Mean pulse profiles and polarization parameters for PSR J1141— 6545 at two 
epochs, 2004 July 07 (MJD 53193) (left) and 2007 August 04 (MJD 54316) (right). In the 
lower panel of each plot the solid line is the total intensity (Stokes J), the dashed line is the 
linearly polarized intensity L = (U 2 + Q 2 ) 1 ^ 2 and the dotted line is Stokes V = I LH — I RH . 
The position angle of the linearly polarized component is shown in the upper panel of each 
plot. Note that postion angles are defined according to the IAU convention; they are absolute 
and apply to the frequency marked on each plot. 
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Fig. 6. — Evolution of the linearly polarized intensity L = (U 2 + Q 2 ) 1 ^ 2 (left) and Stokes V 
(right) profiles for PSR J1141— 6545. The mean MJD for each profile is shown. 
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Fig. 7. — Variations of linear position angle (PA) for PSR J1141— 6545 mean pulse profiles 
as a function of phase and time. Error bars (±lcr) are plotted on every third point. Position 
angles are absolute according to the IAU convention and have been rotated to refer to a 
reference frequency of 1400 MHz assuming a rotation measure of —93.4 rad m -2 (see below). 
For the last few epochs, PAs in the leading outer zone have been offset downward by 180° 
to maintain continuity in the time-variation of PAs in this zone. The mean MJD for each 
profile is listed. 
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Fig. 8. — Evolution of the mean absolute position angle (in the IAU convention) for the 
centre (left) and outer (right) parts of the PSR J1141— 6545 pulse profile at the reference 
frequency of 1400 MHz. The y axes in the two plots have the same scale. 
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Fig. 9. — Time variation of (apparent) rotation measures for the centre (left) and outer 
(right) parts of pulse profile. The dashed line in the righthand plot shows the weighted 
mean RM. 
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Fig. 10. — A representative set of normalised total-intensity profiles for PSR Jl 141— 6545, 
together with the fitted gaussian components at each epoch. The sum of the gaussian 
components is plotted as a thin line, but it is mostly hidden by the observed profiles. 
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Fig. 11. — Time variations of the fitted widths and amplitudes of the six gaussian com- 
ponents. The amplitudes have been converted to component mean flux densities using 
S = * a * w/ P, where a is the relative component amplitude and w is its width, and 
then scaled so that the mean flux density summed over all components equals the smoothed 
value in Figure [2] at each epoch. 
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Observer 



Fig. 12. — Definition of angles used in the analysis of precession of the pulsar spin axis (after 
DT92). The triad (I J K) has K in the line-of-sight direction and the triad (i j k) has k 
along the orbit normal, with I = i defined by the ascending node. Io is toward East and J 
toward North on the sky, making a right-hand triad with K = K. The pulsar spin vector 
s effectively precesses around the orbit normal k at the angular rate Q p . 
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Fig. 13. — Variations of the fitted central PA of the RVM fit (?/>o) f° r the outer parts of the 
pulse profile. Note that the values are defined in the DT92 convention, are absolute and 
have been rotated to refer to infinite frequency assuming a rotation measure of —93.4 rad m -2 . 
They therefore have the opposite sign to the astronomical PA plotted in Figure El The line 
is the expected variation in ip based on the best-fit precessional model. 
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Fig. 14. — Plot of x 2 m the spin-orbit misalignment angle (S) - reference precessional phase 
($o) plane for fits to the outer- zone RVM central PAs shown in Figure [T51 Black and white 
correspond to a reduced x 2 values of 1.0 and 10.0 respectively. Contour lines at the 68%, 
95% and 99% confidence levels around the two possible solutions are marked. 
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Fig. 15. — Plot of x 2 m the spin-orbit misalignment angle (5) - reference precessional phase 
($o) plane around the region of the preferred solution. The greyscale and the solid line 
contours are identical to those in Figure [TH except for the change in scale. The dashed 
contours are at the 68%, 95% and 99% confidence levels for the global RVM fit. The best-fit 
position shown by + corresponds to the variation in ipo shown by the line in Figure [13j 
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Fig. 16. — Final probability distribution for the spin-orbit misalignment angle 5 derived from 
the central PA and global RVM fits. 
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Fig. 17. — Variation of impact parameter (3 as a function time over the observed data span 
based on the final model parameters (Table [5]). The inset shows an extended time range. 
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Fig. 18. — Two-dimensional shape of the emission beam over the traversed region in a planar 
projection centered on the magnetic pole. The inclination of the magnetic axis (a) is taken 
to be 160°, the most probable value from our precessional solution. The coordinate system 
is centered on the magnetic pole and the dashed circles are at 1° intervals in radius. The 
dotted lines show the path traversed by the line of sight with the upper one corresponding to 
the beginning of the data span and the lower one to the minimum impact parameter reached 
around MJD 54000 (Figure fTTl). Since (3 < and ( = a + (3, the line-of-sight trajectory is 
above the magnetic pole. Furthermore, since a > 90°, the line-of-sight trajectory is "outer", 
that is, on the equatorial side of the magnetic pole. The grey scale and contour lines are on 
a logarithmic scale with a factor of approximately 0.07 between successive contours. 
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Fig. 19. — Relative probability density functions for the mass of the exploding star imme- 
diately before the supernova, the magnitude of the resulting kick velocity and the post-SN 
spin-orbit misalignment angle of the pulsar spin (5) after taking into account the constraints 
on 5 from the PA model fitting. The upper row corresponds to a one-dimensional dispersion 
of 50 km s _1 in the assumed distribution of post-SN system velocity and the lower row to a 
dispersion of 100 km s _1 . 



